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SUMMIT 


In the present work, the stability and control deri- 
vatives have been determined from the fli^t test data. Por 
extraction of aircraft parameters measured values of test 
inputs and the resulting output responses are required. Ilhe 
measured responses are conpared with the calculated re^onses 
by assuming initial values of the parameters. She square of 
the difference between the measured and calculated responses 
is minimized by the modified Newton Rsphson method and the 
Davidon-ZLetcher-Powell method. The values of the praaeters 
are updated at each iteration and the process is continued 
till the minima is achieved* The results obtained for a light 
subsonic aircraft by the two methods have been compared. It 
is found that all the important derivatives extracted by 
the two methods, are in close agreement. Olie computed time 
histories by the two methods match well with the fli^t data. 
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aMp m 1 

INMDUGIIM 


1.1 GENMUL ; 

Over the past several years, there had been a renewed 
interest In determining dynamic aircraft parameters such as 
stability and control derivatives, from flight test measure- 
ments'. Ihe need for theee data has long persisted. However, 
only recently hi^ly automated data acquisition systens and 
advanced estimating techniques have been available, that one 
can extract such information efficiently. 

Various methods of estimating aircraft parameters, i.e* 
stability and control derivatives, by use cf wind tunnel mea- 
sur allots or from theory, are in existence. Ihese derivatives 
can be used to conmute time histories of various motions of 
the aircraft. It is important to realize that the conputed 
motions and stability cxe meaningless, unless the equations 
of motion and stability derivatives are truly representative 
of the aircraft -under consideration. YiTind tunnel measurements 
are usually made with small models, and Reynolds number, rough- 
ness, tunnel wall effects, Mach number etc,, generally are 
not properly scaled to simulate the fiaLl aircraft. It, there- 
fore, becomes necessary to assess the correctness cf the dcri- 
vat^Yes, so deteimined by the wind tiainel measuremai ts. One 
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hod is to obtain the derivatives by fli^t testing of 
■craft. In this method from a given groi$) of measarements 
velocities, positions and accelerations taken over a time 
:erval, we have to determine whether the stability and 
itrol derivatives can be estimated. She pro cess of ertr acting 
Qerical values for the aerodynamic stability ard control 
rivatives, is known as aircraft parameter identification. 

The following are the needs for ident ificati on;~ 

lo provide input information to aircraft simulators. 

To provide basis for design of flight control systems. 

The stability and control derivatives define a given 
rcraft nnre uniquely than the response mode criteria (As' 
sted in various flying qualities specifications like MI1“F~ 
85). It is thus anticipated that in future these parameters 
11 ultimately play a major role in the design, testing ard 
rtification of aircraft, 

3 BAQEC-GROIMI) AHl BBIEF HISTORIOAL DEmOPMEKT ; 

One of the first flight test programmes to obtain 
.antitative measurements of aircraft aerodynanic characteris- 
es was reported by lorton in 1919. lift and drag coefficients 
re determined by equating the lift to wei^t aid drag to 
e thrust. 


Soule and Wheatley (e.g%1) ^pear to be the first 
to have determined all the major longitudinal stability and 
control derivatives of an aircraft from fli^t test data. IHtiis 
analysis used siE:plified equations representing one~degree- 
of freedom of motion. Equations were solved for one parameter 
at a time, assuming values for other parameters based on wind- 
tunnel tests. This basic approach was used till middle of 1940' 

Milliken (1947) pointed out that the increasing use of 
automatic control systems required more accurate modelling of 
the aircraft dynamic characteristics. These factors coupled 
with the research engineer's motivation to inprovetbe accuracy 
of fli^t results, stimulated the development of several new 
techniques for determining stability andcontrol derivatives 
from fli^t data. 

In the late 1940* s throu^ mid I950's servomechanism 
theory was expanded r^idly and the frequency -domain techniques 
of Hyquist and Bode (see e.g.2>3) were popular. A disadvantage 
of this approach was the considerable flight time required to, 
sweep through all the frequencies of interest at each flight 
condition. 

The problem of determining stability and control deri- 
vatives is based on a linearized, small-perturbction mod^ of 
the aircraft dynamics. Hence it was natural to consider using 
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a linear least-square fit of flight data to the linearized 
equations of notion as was done by &’eehhi9£g(4) in 1951. In 
1954 Shinbrot (5) developed a geaeralized least-square nethod 
which encongjassed the earlier least-square methods* During the 
period when Greeuberg and Shinbrot developed their DEbhods, 
the digital computers were not available, fhese two methods 
required extensive calculations to match the computed response 
with the flight measured response of the aircraft, Ihus a large 
volume of flight data had to be processed manually for extrac- 
ting the stability and control derivatives. 

‘inalog Matching’ (see e*g. 6,7) technique was used as 
early as 1951 to check aircraft parameters determined by other 
methods* ©ren though the technique of analog matching has grea- 
tly ic5)roved over the period, the accuracy of the ©suits is 
highly dependent on the skill of an individual operator* If 
several parameters are to be determined, this technique requires 
excessive number of man-hours to obtain an acceptable solution* 

Jlthough several attempts were made throu^out the l950’s 
and early l960’s to inprove techniques, the effort was rela- 
tively small. Two factors caused a revolution in aircraft para- 
meter estimation techniques starting in Eiid-l960*s. These weres 

i 

1. Highly automated data acquisition systems were becoming 

standard in flight testing* 


2. lorge-capacity, hi^-speed digital computears were 

available to solve cocplicated algoritbins efficierrOy, 

The ability to transfer the flight data directly to the 
coEtputer with no nanual operations on the data, helped 
in increasing the Use of fli^t-testing techniques for 
evaluating aircraft parameters. 

The interest in parameter estime-tion was renewed in 
1968 . iarson (8) applied the method of quasi-lin eari2 ation at 
The Cornell Aeronautical Laboratory, Taylor ^d Iliff ( 9 ) used 
basically the same ^proach, bub referred to it as the modified 
Newton-Raphson technique, at the RASA Rii^t Research Centre. 
Denery (10) applied the general theory of system identification 
to air vehicles for extracting aircraft parametera,Melira (11) 
used the method of stochastic approximations to solve gysten 
identification problems, later on the same ^proachwas ^plied 
for aircraft identification (see 13). 


1.4 AlRGRi.RT PiUilMETER ESTIMiiTlOR ; 

The general problem of parameter estimation is to 
determine certain characteristics of the physical system from 
e3q) eriment al test data. Measur^ents are made of the test in~ 
puts ^'d the resulting output re^onses that depend in some 
way on the systen characteristics to be determined. Paraaeter 
estimation technique is the process of estimating diar act eristic 

I 

from tlie trput/output measuranents* 
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* The general concept of the aircraft parameter esti- 

mation technique is illustrated in Pig. 1. The various aspects 
of the problem are: 

(i) Mathematical Model 

(ii) Estimation Criterion 

( iii) Oocp ut at i on al . gorithia- 

(iv) . Instrumentation and Da«.a inquisition 
(v) Test Ii^uts 

1 • 4* t Mathaa.atical Model ; 

A mQ,.del must be selected that adequately represents 
the -aircraft <^aract eristics to be- measured*' ‘Tbr the -present 
problem the aircraft is assumed to be rigid and the linearized 
equations of motion are considered. However, in other instances 
a more coEplex model may he necessary, such as. at very high 
angles of attack a non-linear model may be required, in inappro- 
priate- model can degrade the accuracy of the parameterrbstimat es. 

1*4.2 Bst imati on Criterion ; 

There must be some means of assessing fit of the 
conputed re^onse to the measured response. To irpleaent this 
* criterion function* is used. It is usually some form of the 
integral square of the error between conputed and measured 
response. SJhe best estimate of the paraaeters is the set of 
parameters that minimizes the criterion function. 
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1.4*3 Qomputational .dlgorithm ; 

The criterion function is often non-linear with 
respect to the parameters to be estimated; therefore it has 
to be minimized by an iterative computational algorithm. 
Important factors in selecting the minimization algorithm 
are convergence, computation efficiency, local minima etc. 

1 . 4 • 4* Instrumentation and Data .acquisition System ; 

Aircraft parameter estimation is highly dependent on 

the quality of the fli^t measured data. In reality bias and 
random errors 'can arise from improper location or orientation 
of sensors, calibration of measurement and recording system. 
Other errors can be introduced from electrical noise, engine 
vibrations, inappropriate signal filters, air turbulence etc* 
Any elimination of errors, noise or 'unceri; ainties within the 
data acquisition system will ic5)rove the accuracy of the 
estimates. A comprehensive discussion of flight-test instru- 
mentation for aircraft parameter estimation is available in 
Ref. 12. 

1.4.5 T est Inp ut : 

As a minimum requirement, the test ii^jut must excite 
the principal response modes that depend on the parameters to 
be determined. In the example considered in this work, a combi 
nation of rudder and aileron piCLses adequately excited the 



8 


later al-airection raotion. However, a q.uestion arises 

wH ether one t^o of control irput ni^t be better than the 
other, in the sense that it provides better estiaates. 

Several ps^ers have considered this question (15) but the 
concept has not been fully explored in a flight -test appli- 
cation (14). 

1 .4 PRESHM'f wore: ; 

In the present v/ork, we are interested in analysing 
the flight data, in order to obtain the stability and control 
derivatives. Ihe details of flight data acquisition and those 
of test ii:5)uts have not been considered here* A. subsonic 
light aircraft is selected for the analysis in the lateral- 
directional oode, the flight data for the same coiad be 
obtained from Ref. 17* 

Ihe process of extraction of stability and control 
derivatives involves the measurement of test inputs and the 
resulting output responses* Ihe measured responses are com- 
pared with the calculated responsesby assuming initial values 
of the parameters. Ihe square of the difference between the 
measured and calculated responses is minimized by the modified 
Newton Raphson method and the Bavidon-iletcher-Powell method* 
Ihe values of the parameters are updated at each iteration 
and the process is continued till the minima is achieved. 
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IThere is no direct referonco in the literature about 
the use of the Davidon-i'lctchcr-Powcll aethod for the flight 
dynanics use and its results. It was, therefore, decided to 
use this method to find out if it could he f^plied for solving 
aircraft dynamics prohlesus. 

Ihe details of .the methods for determination of 
stability and cn ntrol derivatives from fli^t data are pre- 
sented in Chapter 2. Ihe identificotion problem is stated and 
a brief review of the various methods used for aircraft para- 
meter identification has beeri carried cut. The modified Newton 
Raphs on method and the Davi don -Hotcher-Po well minimization 
techniques are discussed in detail. 

The corput ational details, results, discussions and 
conclusions are given in Chapter 3. The parameters estimated 
by the two methods and the time histories are compared. Finally 
the source of errors are pointed out and the saiggostions 
are made for further work. 
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GIUPIER 2 

MEIHODS gOR MMIMIOE OF STJ^BILITY M 
nolKiOL DmiVuJ^IYBS IhOM FJjiaET E-k TA 

2.1 IITTROEUOIIOH: 

A general statement of the pro 11 era of aircraft parameter 
identification, as ^plicd to the extraction of lateral~dir ac- 
tional stability and control deriTatives, is presented here. 

A brief review of the previous methods is presented first. 
Subsequently, the two techniques chosen for the present 
study viz. the niodified Hewton R^^hson technique and the 
Davidon-Eletcher-Powell minimization techniques are described 
in details. 


2 , 2 SliiSlLllY lElEiMBIiJllOH; 

STAUmmT OF ffliE PROBLEM 

Por most of the aircraft dynaaics problons, the longi- 
tudinal mode is usually seperated from the 1 at eral -direct ion 
mode, because the resulting error is small. She estimation 
of stability and control derivatives in lateral -direct ion 
mode is more complex than the longitudinal mode. Shis work 
deals with the problem of stability and control derivative 
estimation in the lateral-direction mode. 


!I!he mathematical model often used to describe lateral- 
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directional 

^airplane dynamics can Ic expressed as a sjstesj of linear 
ordinary dif‘f erential equations v/ith const ant -coefficient b> 
in tile following form (see ij^pendix j*) 


Jc(t) *= iix(t) + Bu(t) 


where 


x(t) 


r .1 


I*-* 0m 



p i 
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r-' 0H 


f 
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4 

x(t) = 

r 

II 

a 
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1 
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( 2 . 1 ) 
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0 
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i 

0 

0 
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3h.e last column of matrix B represents the effect of uncertain 
Bias on the measuroaent of i), f and p • 

let z(t) and 2(t) be the noise contaminated measure^ 
ments of x(t) and xCt) respectively* !Ehe control ir^uts 
u(t) are considered to be noise-free* 



n 


Given z, z and u, !:iie problen is to deternine eertain 
■unknown elenents of natrices A and B, i<,e. the stability 
and control derivatives. 


2-5 BliEVIQUSJiMl iO]^ -, i .MLBSJiMLm i 

Some of the nethods used earlier for deteroining 
st’ability and control derivatives are described below: 

2*3 *1. Simplified Equations Method : 


. for selected type of responses, the effect of only 
a few coefficients dominates, thus permitting thedse of 
single e^qpressions to deternine these coefficients, lolowicz 
(12) obtained good approximations for some of the longitudi- 
nal stability derivatives by keeping only the dominant terns 
when equations of motion have been solved for a particular 
derivative, for example, G can be determined from the 

jlJ o 

^ e 

initial portion of rapid pulse maneuver by s 



( 2 , 2 ) 


Similarly, 
from the 


G can be approximated to within 5y. accuracy 
“a 

relation. 





1 i 

2 



1 

2 

V SO 



s.p# 


(2.3) 
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where W is the frcouenoy of short pe-riod oscillatiouc 

“s.-p 

Wolowicz (12) c-lso £:ave some simplified ^proxinat.. , .is 
for the lateral stability derivatives. However, because 
the more complex behaviour of the airplane and large nun!. ■ 
of derivatives involved, the lat eral~dj.rectional control 
and stability derivatives are not as readily and reliably 
determined by the use of s^proximate equations, as are thi 
longitudinal derivatives. 

The disadvantages of this method are: 

1, Only some of the primary unknown coefficients of stab-' '*7 
and control derivatives can be determined, 

2, The forms of response that can be analyzed are very 
restrictive i,b, effect of control must either be dominant 
or negligiblet 

2 .3 V 2 * inalog-'Matchin^' Method : 

In analog matching technique, the values of stability 
and control derivatives are assumed and motions are computed 
for the same control inputs as tliose used in fli^t tests^ 
Comparisons are made between the computed aid flight-meas" '-d 
responses. If ary of the responses do not compare favouroiy, 
some of the derivatives are changed and the motions are 
coB^juted again. The process is repeated until satisfactory 



agreeingnt is obtained between the estiiaated and flight - 
measured responses, fho ability to converge to a set of 
acceptable derivatives, or of even finding acceptable agree- 
ment between coiitputed and fliehb notions, is largely a natter 
of experience. 

The main disadvantages of this technique are: 

N % 

1. ibnalog-matching depends quite heavily on the experience 
of the operator. However, this method can be us ed as a 
back-up method and has given satisfactory results in 
many cases. 

2. The method works successfully only when a single control 
surface is moved at a time so that the maneuvers are 
simple, (see 12). 

3* When the maneuvers are made with stability augmentation syst 
em, the data is difficult to analyze. 

4. This method is extremely time consuming. 

2.3.3 Least-Squares Metho d: 

Least-squares method assumes a performance criterion 
which minimizes the square of the state equation error, by 
substituting the measured value of the state and its deriva- 
tives i, e. it minimises the following cost functionals 

J = / (2; Az - Bu)^ (z - Az - Bu) dt 

o 


(2.4) 
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To derive expressions for the values of ^ and B that 
nininizes Eq. (^.4), let 


0 = C ^ i 3 -J 


Then, 


T 


T 


^ :p J Z* Z dt 
O 


T 

2/2 

o 


[ ] dt + 

U ^ J 


1 z 

/ r- ~ -1. 

; u i 


T ' r Z 
c c [- 


u 


-] at 


(2.5) 


The niiniiiiizafcion is achieved hy talcing the first derivative 
of Eq, (2.5) w.r.t. c i, c. 


This yields, 
_T 


9 c 



-] L“ 


Z 

u 


T 

-1 * 


-1 


j 


T 


/ [ 1 

; L u J 


•T ,, 
z dt 


( 2 , 6 ) 


Eq, (2.6) is the desired solution by the method of least- 
squares. It has the advantage of compact form hut disguises 
the independence of qach of the equations to lie minimized. 
This can be shown in the following manner, Consider only 
the first state equation from (2,1), 

«.j(t) S= x(t) + u(t) 
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where 

\ Ip 0] 

Si = 

t 

Gi = C j 3 


Then 


I 


J = / Ti* 


r-, ~ L P 

o "^e P 


- JLi_ r 
e e 


\ B - L, .6 


L. .6 - L 1 dt 

r oh 



where the subscript denotes the neasured value. 

e 

Once again taking the first derivative and setting the 
resulting expression to zero, gives 





( 2 . 7 ) 


Now G. is the first row of the matrix c and p„ is the 

• . 

first eLement of z, which makes it apparent that the elements 

of the first row of the c matrix are indep enden-^f all the 

elements of z except the first, p^^^ . it similar relationship 

e 
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can be 


easily shown for th’- 
The disadvantacos 


of c matrix. 


1. The row independence (ohcwn vc Ly Sq. 2 . 7 ) is one 
of the drawbacks of this ncxh-,*!* in Ihr.t only one of the 
measiired state variables is us*...: in dotcrnining a given row 
of c matrix. If one of the eigh-J-o has not been measured, 
the least-squares method dots n.t i^rovide ai estimate of 
the derivatives rd.ated tc that sirhal. 

This independence also illuctrntcs that the estimate of 
one row of the c matrix is obtruncd in.lcpcndently of the 
other rows, and no *trade-off * cmri tc made between elements 
in different rows to improve the match. 


2. If two or more of the mcasur<-.; responses are linearly 
rdLated, this method gives an ahsuri; ' solution or no solution 
(15). 

3. This method gives excessive variance of the estimated 
coefficients (16). 


2.4 MOilPlED NEWTOIif-IUPHSOI MVIPON-gLETQHai-POWELL 

techniques ; 

2.4.1 General i 

These techniques are ]QQ«nii a£ sheeting those para- 
meter values whida best fit an assumed model to a data set 
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according to a particulor error criterion. The error criterion 
is more general than the least-sciuares criterion in that 
it permits the fit error to p, r, p and ^ to be minimised 
as well as the fit error to p, r and p. These techniques 
also enable one to use a priori values of the stability deri- 
vatives, bias terms, and initial conditions to in^jrove the 

match 

fit of the equations toZfli^ht test data. 

The mathonatical model chosen here describes the lateral- 
directional motions of an aircraft by linear’, constant coeffi- 
cient differential equations. Consider the following model 
(see Appendix A) 

X f= Ax -t Bu (2,8) 

and y = [_ I _] X + C- g ^ u (2.9) 

where, 

The vector y is the set of out-put response quantitiesj 
y=[;prp(j)prp*] r 

The matrices -i*-' and B are defined in Eq, (2.1) 

' [ I ] = Identity Matrix 
0 J - Null matrix 
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Let 2 denote the measnremenit of the actual aircraft 
response quantities corresponding to the computed output 
response quantities, y* The measured response s woiiLd not 
he exactly the same as y because of, 

(i) The measurement errors 

(ii) The difference between the actual and the assumed 
linear mathematical model, Eq. (2.9) • 

3he objective is to minimize the difference between z and 
The appropriate cost functional is, 

T 

J *= / (z-y)^ D-,(z-y) dt (2*1C 

o ‘ 

Where is a weighting matrix reflecting the relative 
confidence in the measuremoats# 
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specifying the cost functional J, an algoritlin 
is to he chosen for minimizing this cost functional. Many 
methods are availahlo for non-linear minimization, Hov/every 
for this particul.ar work a modified Mewton-Raphson technique 
due to Hiff and laylor (9) and the Ravidon Fletcher Powell 
technique are selected, Ihey are described below, 

2.4*2 Modified Re^ t on-Raphson Minimization technique 

For convenience, let us define a column vector c 
of the unknowns to be estimated. Ehe elements of c are some 
or all of the unknown elements of A and B of the noise biases g 
and of the initial conditions XjCo), e.g. 

In all subsequQit calculations, the noise biases g^^ and the 
variable initial condition (o) have been neglected*. 

The Newton-Raphson technique is an iterative method 
for finding a zero of a nonlinear function of severe para- 
meters, or in this instance, a zero of the gradient of the 
cost functional J, i, e. 

Consider a two -term Taylor *s series expansion of 
Vq ^ about the k^^ value of 
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(V„ (VeJ)k+ (V%)fcAo^, (2.11) 

wjaere; 

'^^‘tc+i “ 'Vi ■ 

tils scGontj. gra.dicsi'fc of "tliG cost functional 

w.r.t, at the k ^ iteration. If equation (2,11) ig a 

gy.fficiently close approxifiation, the change in c on the 
"t/ ii '' 

(k+l) iteration to nake fV^ J)] 5:^1 Ejjproximat ely zero is, 

AVi=-C(VcJ)k3 (2.12) 

Ihis i#” ITewton-Saphson algorithm, 

lo evaluate the chrhge in parameter value A c^^ by 
Eq. (2,12), thq Voalues of (% J)^ and are needed. 

She gradient of J with respect to the vector c can be expres- 
sed in 1y«ifms of the gradient of (z-y) by dif f erenti ating 
Eq. (2.10) as, 

C2,13) 

In order to evaluate ^ the term SJ^iz-y) ±s 
required to be specified in addition to the terms, already 
defined. Ihis can be obtained by differentiating Eq.(2.9), 



ns iQllcmut 


U “ - (Vo C- 5 -]V - C- ^ - {^0 C- i -]} 


- C- 5 -iVc 


U 


( 2 * 14 ) 


(z-y) om be expressed in terms ojT various 
partiiii derivatives* 33aese partiai derivatives with, respect 
to the individual coefficients of c (^i|» 


I 

M- s -I _ 

• o o 

I 0 , 

81-5-1 __ ■ 

0 . 
0 

• 
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* < 

0 

mm 

I ■ 

•l-5-i _l 

0 

0 

a j- g “{ 

0 

• 

9 

8 

• 

0 

p.. ,, -[X = 

8 

T 

% 


where x^ aod u^ appear in ii-p row of the vector^ 
The isradient of x with reference to c is given by 
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Vo* = 


9P 

ap 

ap 

a c.| 

tC 2 

a c. 

ar 

9 r 

ar 

• • • 1 -,- 

a c^ 

8 Cp 

ac. 


..li. 

ap 

t # # # T r 

a 

ac2 

ac. 

dp 

a<I) 

. . . . ±i 

a c^ 

ac2 

ac 


n 


n 


n 




The elements of grrflent x Cc?n be determined in the 

"t bi. 

following manner. Expanding x at (k+l) iteration by laylor 

Ic 

series, about x and retaining terms upto first derivative 


,.k+1 


X 


k •k . 

X + X ^ 


m 

Substituting the va^uo of x from Eq_. (2.1), we get, 

x^^ = xVaC ^ ] C'2.15) 

Taking partiaL derivative cf Eq.. (2.15) with respect to 


a^^ and i't becomos, 


r^x 


k+1 


M 


k 


®i3 


6 a,. 






aC-^ 


''=13 ®=13 


k+1 

fix 




ab 


13 


'^3 


ab 


13 




( 2 . 16 ) 


The initial conditions are invariant with respect to a^ • and 

ij 

bjy, i. e. 



2 -:- 


1 

a 


2£l2.1 



0 


ancl 


0 


'Tile elonents of the gracliiait x ca:n now be determined by 
Eq. (2.1 6), im alternate procedure to determine the elements 
of gradient x, due to laylor (9) is given in Appendix B. 

T.hus all the terns in Eq. (2,14) btave been evaluated; 
XJ J can now be evaluated by Eq. (2*13), 

Ihe difference between the measured and confuted 
responses (z-y) , can be represented as quasi -linear with 
respect to a change in the unknown coefficients, i,e. 


C 2^-y 1 k r 2-y ] ic_i + Vo [2-ylk^“k 


Using this approximation in the cost functional results in 
the following first and second gradientss 

Vc C Cz-y)j^ ] (z-y) 2k dt j- (2.17) 

y/j - 2 / (z-y) CVe (2--y)3k dt + 

1 2 ® 

(z-y)^[ (z-y)k dt (2*18) 

2 

Ihe second term ofv^ Jdininishes, as the response error 

o 

(z - jj^) diminishes. Ihe laodified Hewton-Baphson technique as 


;■) 


SSI 


t tir'.t. 


'ii n 


R' 


25 


developed bj Taylor (9) ne{.''lc-cts this tern. Thus, the final 

2 

relation for^ J hecones 
o 


2 

Vo " 


C (z-y) ]j, dt (2.19) 


Now, substituting the values of "SJ? J and V J ^1* (2.12), 

' c o 

the nodified Newton-ilaphson algor itlin bee ones 

^“k+l " '{o “'f ^ * 

/{[Vc (z-y) Jic} (2.20) 


ill the terns in Eq, (2.20) involve only the first 
gradient of (z-y) . This greatly reduces the computation 
time and the approximation intproves as the solution is 
’approachedf 

For computational purposes, the integrals are treated 
as summations. In the indicial notations the Bq. (2*20) then 
becomes. 

Ad-- CVo %Vo V) } * 

3-1 


( 2 . 21 ) 
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wliere eul3sc]pipts % end j c.re the indices indid^ting the tine 
sample, end 1 is the total humhon of sacipies* 

Equation (2,21) is the final form of the modified 
Newtdn**Iiaphson algorithm 


2 . 4 • 5 • Day i do n-ZL etch or "-Powell (3^»E«P«) Minimigation lechniaue i 

2 ' 

In the D.F.P. technique the local aeesian matrix (Vc 

Ic 

is replaced by an qjproximete matrix She method of conpu- 

f' 

ting this matrix co;iq)letGly ^iminates the need for evaluating 
second derivatives ani performing matrix inversions, and yet 
the sequence of iterations converges quadratically to the 
minimum point (18)* 


ifter selecting a starting point, a direction of search 
is coEiputed as follows (i9)? 


Mi 



( 2 . 22 ) 


where i, ;| , 1=1, 2 .... m, m being the number of unknowns 

c., and k is the iteration index. M are the direction vector 
^ i 

ft J 

coiaponents, gradient vector components and H. 

D % 

are the elements of a positive definite matrix (mxm), which 

is initially chosen to be identity matrix^ It is evidQ^t 

from Eq. (2^22) that the initial direction of the search is 
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the path of steepest dcsccKt, 

A one-diaensional S’-arch is conducted in the direction 
chosen hy Eq. (2,22) until c nininuio is located utilizing 
the relation, 


k+1 

c 

i 



* V 
a 


( 2 . 23 ) 


where a is the stop size- in the direction of search. 

how a convecgenco check is nade. If the convergence is 
achieved, the process is terminated, otherwise, the matrix 

Ir , -I 

H is calculated as follows; 


1?+^ = 


where 


Ac^ ( Act 

(Ao't® (( AS)*") 


T 


(AG)^ r(^Cr)^1 ^ 

(AG)^ (AG)^ 


k k+1 k 

A c c ~ c 


( AS)*"= 


(2.24) 


Ihe value of updated matrix H is substituted in 
Eq. (2.22) and a new one- dimensional search is made. ®ie 
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pro cess 


is repent ci till the conTc-rgence is achieved, 

. (2.22) to (2.2-v) is the D,?,!*, algorithn (see19). 





■D.BI-iJ I LS. RES UL TS iilD DISCUSSION 

3 . 1 Iii£iSQJ^£IOl3j 

The fli£jit test (Irt.- for the Ir.tcral-dlrect ioncl loode 
of a li^t subsonic aircraft, flying at a speed of approxiuat ely 
280 M.P.H. , arc obt-oincl fren Ref* 17, i'his data is used for 
extraction of 1 at oral- direct ional stability and control deri- 
Yatives by the ixjdified IJewton Raphson method and the favidon- 
Zletcher-Powell method. iSie computational details are discussed 
first. The stability and control deriYatives obtained by the 
two methods are coraporod. The aircraft time histories calcu- 
lated by using the extr nctj cil derivatives are compared with the 
flight test data. The results are then discussed and the sources 
of errors are listed* Finally the conclusions are drawn and 
some suggestions are mace- for the further work. 


3 . 2 GQBCPUTid^IONlL DETJhS : 

The flow diagram depicting the steps involved in the 
extraction of lat erol -direct ional stability and control deri- 
vatives from flight data, is given in Figure 2. The listing 
of the programmes developed for the modified Newton Raphson 
method and for the D.P.P, method ore given in appendix C and D 
respectively. The different variables used for the computations 
are defined in the program listing* 
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I’ii.., iri,..ut t 'r'f ,„rc^j,nus cciisist of the initicl 
^niesscu .-f thu vit oio contr':.l 'derivatives, the control 

input state, n.'i th'.. •or :, it .Ou.-Tsai'o:. responses p, r, p, <f), 

* 

p, r nnl p. Ik r;boaio^ .•i‘ ’liie specific inf'ornation onl for 
convenience, th^ v;..iijatin,' •■.r.trix is tchen '-s identity 
natrix* Since the volacs of y^ and y^ are generally snail, 

3? 

they are taken as zcT'.. (sec e. g. 16,21). Ihe value of y^j^ is 
calculated fron the speed ,jf the aircraft. Ihis reduces the 
nuniber cf unimowns t;.; 15 i.c. 8 in natrix A and 7 in natrix B. 
Ihe responses ere neasured frtn 0 to 10 seconds at interval 
of 0,1 seconds. Eius there are 101 data points of aircraft 
responses. It is asounod that only the values of p, r, p and '■ 
^ are known. The values of p, f, and ^ are calculated by 
differentiating p, r and p by finite difference nethod. Ihe 
tine histories input is 'piver in .^jpendix B. 

She Newton Raphson iDrcgrurne consists of the nain 

programe and subroutines Giilhli’ and M^tJEINV. The subroutine 

GhLEEP calculates the rosrponse of the aircraft and (z~y)» 

denoted by OIK and H in the progracne. The subroutine M1IIN7 

2 

calculates the inverse of \7c J. The increnent in the deriva- 
tive values is calculated in the nain prograa and the values 
are updated at each iteration. The process is repeated till 
the nininun value of the cost functional J is achieved. The 
final values of aircraft responses are given as output in 
■iippaindix 
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fhe 13,1 .P, propgraine ccnsists of t}je main programme 
$n<i subroutines Q^JKf, PUKCT and DSP, Ihe subroutine CULDES 
calculates tho response of the aircraft and y ^ # Ibe 

subroutine SUIIGI calculates the value of cost functional J 
and its gradient VeJ. Ihs subroutine DiP first carries out 
a one dinensionol search until a niniiaa is located. If 1he 
convergence is achieved, the process is terminated, othervase 
the value of the matrix H is updated and a new one dimensional 
search is made, Hu. values of theparamet ers are updated at 
each iteration, Hio iterative process is continued till the 
DiniEia is achieved. Hic final values of the aircraft responses 
are given as output in j^ppendix G-. 

3*3 RESUIfS MD msCUSSIONg 

The aircraft time histories calculated from the stability 
and control derivatives, extracted from the flight data, by 
the modified Newton R^hson method, are shown i'n Pig,3> 
time histories obtained by the D.P,!. method are shown in f'ig,4 
It can be seen from these figures that in both the cases the 
aircraft time histories match the flight data well* The fit 
obtained by the method is slightly better than that 

obtained by the modified Newton Raphson method. It can ^so be 
noted that none of the responses Mve ^Own any difeb^^ce 
from the fli^it data* There is a slight viatiatton at the pa^ 
values of some of the responses# Eie possible reasons the 
same .are enweraled in Section 3*4’ 


The, ■/' 


!:t Cv-ility rjid control derivatives 


. 0 , lisred in Table 1. It is 
tnat v;d.u ei;.,ht .)i the twelve derivrti'VGs 


cbtrdn.,d by the. tv:, 
observed 

extracte,*: fr i:: the o‘li •;.t ly the tv/o nethods, ere close 

to each other. Thu vru-u.-s h' tiic * dri^.-in^; -in-roll* derivative, 

L ; the siJx force .leriv' i: c dUu to side-slip, , and the 

P p 

yaw due t roll /lorivativu, IT axe sufficiently in error. 

The sigh of the roll to yf.-;/ derivative, is reversed. 
One of the reasons of this co-uld bo that these parameters are 
v.:ealc parrmetex:^ that cannot lo accurately identified from 


fl ight t e st d at r', (15)* 

It has boon itliowri by itcderick (20) that stability 
derivatives I,. > ^ j II ni'c m;'. st influaitial in deterinining 

P p 

lateral-dirocti .nal. stability. The table 1 shows that 'tiese 

derivatives have been recovered accurately by both the methods. 

It is also noted that even though some of the unimportant 

derivatives such as IL, Ioj sufficiently in error, the 

p p r 

time histories rxc closely matched with the flight data in 
both the cases. 


The modified Newton Rophson technique when used with 
zero initial values as starting values of the. stability and 
control derivatives, gave large fit errors* in erroneous set 
of stability ^d control derivatives was obtained and the cal- 
culated responses differed considerably from the measured 
aircraft dynamics. It is, therefore, essential to provide 



a^'T’-pri -vtT^ ir;i*.iai 

c -■ *' ■, - **,ii 4,. 


‘ ‘-icr; 


■t . t:icsc JcrivativGs wldLch. 
i- C-iic r.iatric'wS x.. and B, 




'lan: to be sensitive to 


tile initial ..;t ".rtin, .•.-i 
derivatives, 

i’hu .o,;::'ii i,.; . i;.,',;- 
to n rdniinu..-. vrduc .!,* V... 
the nah,;. never r, 

alv; ry 3 c o nv er <; -s. t i 
can be satisf aot .:.ril v ur;, 
Raphson techni.jue I 'die i 
of the iholific: Id-.vtjn u- 
nethod, resultiin inl.,;:o 


• ‘ .. till, stability and control 

iiaphe.sn teclmique first converged 
.■'i’i error and then diverged. However, 

.. ^.e.: this tendency, fho D.B.P. method 
■ml r-lnima. fiius tho D.P.P. method 
f’- cases where the modified lewton 
aonvurge. fhc convergence in case 
ih'on muthsci ivas faster than the D.P.P, 
'4 . u.t vatxonovl t ime « 


Oolculaticnn i’-r tht 


problem considered are performed 


on mfk 7044 


K.V syctc. 


the order of magnitude of the 


computer t.imc rcQuircd l,y tjit. two tochniq_ues for solving the 
present problem, is as fullows: 


Modified licwton idphsen technique 10—12 minutes 
Pavi don-PL et chcr-Pov/ ell technique 3 5'"40 minut es 


3-4 SOURCE OP 


She aircraft behaviour cannot be predicted very 
accurately by the linear dynamical model of the aircraft. 
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<'2id ilu: 


H' c:.'u 


'uplin^^- between the longi- 
^arvi-uiri.ciieni.xL is accurate 


en.'iu^i for tiu fiirr.t 
this the fact- re t::r 
pute,'. rfflu r.i., r’iiure ti; 


':f the nnalysef clata. Besides, 
fjr the difference in the con- 


iuatcricvS 


cf the aircraft 


apd in the 


evaluation ..f otaiiiity 


ccnurjl derivatives, are as followsi 


1, I12ie inj;. ut piiis,.: r. 

in the viclr.ti..'n ,f a;ofu:.: 
sing r-io^slitude will i n cr- 


ib lit u^c is not kept saall resulting 
Ji-ticns cf saall perturbations, Increa- 
.uc c inertial :r aerodynauie non- 


linear it ii s. 


2. Only tnc p'U'ti'.il ac.tr v/as av.'iilnble for the analysis, 
Ihe value .,:f the reeponroj variables p, r and p were calculated 
by differentiating p, r and p by the finite .lifference method. 
Ibis intro tiuces error in tiu. estimation of p, r and p and thus 
in the estin.ati ni of paraivictcrs and the aircraft tine histories. 

3, Since the ri. suits cf all parameter identification 
procedures depend heavily on the quality of the test data 
avail £l>le, there is a need fer min ini zing the instrumentation 
errors. Some of the mrc conmon instrument induced errors 
include random noise, calibration errors, mounting in accura- 
cies, instrument bias and tine lag. Consequently, to achieve 
reliable results, the data must be conditioned by coripensating 
for instrument shortcomings. 
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L-: Cl§_ , lie 
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pcroxicxors like vary sufficiently 


no t p r clorfln an t . 


3.5 Gouaiasions : 

Both the tochni;jUcs via. the oodifiod Newton Eaphson 
method and the B.F.P. method £lve good results for the example 
considered for a lic^it subsonic aircraft. Thou^ the values 
of unimportant pcirameters like N^, vary sufficiently, 

their effect on the time histories is not predominant. Hie 
calculated time histories of the aircraft match well with the 
flight data for both the methods* 

Both the methods axe found to be sensitive to the 
initial values of the derivatives. 

Ihough the modified Newton Baphson method converges 
faster than the D.S'.P. method, it i^ows a divergence after the 
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nitii-.d iij riwiiit,\) c,.,* "tiie nciJliocl ttsver shows 

such behaviour* It fa-.vo.ys cclivargoo steadily to the nlhinun 
value* Thus thi:; rjcthoC ccuid prove suporiar in cases 

where the i'iyvrtd li keijhson nothod does ni/t show convergence* 


3.6 

1* /ji algorithm bo evolved which would be insensitive 
to the initial ostinetes of the parameters and con- 
verges fast to the correct values of the perameters. 

2. 'ihe nodified Newton Eaphson method and the D.F.P. 
method should be used to extract the stability and 
control derivatives of high perfoniance aircraft and 
the wide bodied transport aircraft, to confirm that 
those methods could be applied to common conf igur at ion; 
of aircraft. 

3* ill effort be na.de to find out a suitable weighting 
matrix to avoid divergence after converging to the 
minima in case of the Ek)dified Newton ik^hson method* 

OEhe effect of neglecting the second term in Eq. 
(2*18) lie explored to find out its effect on the 
overall results of the modified Newton Bu^hson 
technique. 

ij 2 effort he made to modify the D.E.P. method 
suitably to reduce the cos^ujfcertiion time* 


4 
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i. 


04 MOriOH t 


The lateral -directional equations of motion of an airplane 
.'’xo (21); 


-ip p 

4 p 

■*■ p — Y" " ' r — 1 r — lie ^-r» 

P ixx r ^ ^ 

(j-1) 

-^’p p 

^ZZ 

p - IL p + i - Ij, r = 6a + KjBt 

3i 3? 

(A.2) 

p - l( 

3 ^ ■ 

P ijr P 6^ 3- 

(A.3) 


Neglecting the product of inertia Ixz and assuming = 0, 
these equations reduce to the followingi 




(A.4) 


r » p + ^ P + ^6^ 

p - Y p - r + I„ p + \ Sg, + ^5 


where ^ j ^*3. I<q> are the effect of uncertain 


bias on the measui'oaent of p, r and p. 


Realizing that (|) = p . Eq. (■4«4 - ^t) vdaen written in the 


matrix form become 



Equation (i.7) can be represented as 

x(t) *= L x(t) + B u(t) (A.j 

vdiere 



Hote that x, x and u are tine varying while the elaneirts of 
natrices A and B are constant. 
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i^PPMPlX B 

DBPZR.>II:r.4llOIv CP ELEAMJS Qg x: 

■Sie elcccintB of gradient x can be determined in tlie 
following manner. 


aij 


Biffercntigting state Eci* (2,8) with respect to 
viQ obtain 



By solving the differential equation (B*1), we obtain for 
8x 

93J. 


iigain for coefficients in B, we obtain 



B{t- X) 

e 



uCx ) dx 


Siniilarly if the change with respect to initial conditions 
Xj_(o) is considered, we obtain 


dx „ „-^t 
(o) 


G 

0 

6 


, . .th 

where 1 appears m i row* 
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